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Abstract. The relativistic kinetic equations describing time evolution and space dependence of the density matrices 
of polarized photons and electrons interacting via Compton scattering are deduced from the quantum Liouville 
QQ . equation. The induced scattering and exclusion principle are taken into account. The Bogoliubov method is used 

' in the frame of quantum electrodynamics. The equation for polarized radiation scattered by unpolarized electrons 

is considered as a particular case and is reformulated in terms of the Stokes parameters. The expressions for the 
scattering amplitudes and cross-sections are derived simultaneously. 
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1. Introduction 



Compton scattering is an important physical process in many astrophysical systems, such as active galactic nuclei, 
X-ray binaries, and pulsar magnetospheres (see e.g. reviews by Blandford & Scharlemann 1975; Pozdnyakov, Sobol', 
& Sunyaev 1983; Poutanen 1998). Inverse Compton scattering (i.e. scattering of soft photons by hot electron gas) is 
believed to be the main mechanism of the X-ray/gamma-ray production in accreting X-ray binaries (e.g. Sunyaev & 
' Titarchuk 1980; Poutanen & Svensson 1996), gamma-ray bright active galactic nuclei (e.g. Sikora et al. 1997), and 
possibly gamma-ray bursts (Stern 1999; Ghisellini & Celotti 1999). The hard X-rays can interact with the cold material 
via classical Compton scattering losing their energy and producing a cutoff in the "Compton reflected" spectrum (e.g. 
George & Fabian 1991; Poutanen, Nagendra, & Svensson 1996). Induced (stimulated) Compton scattering becomes 
very important when the radiation brightness temperature is large fceTb/mc^ ^ where tt is Thomson optical 

thickness. Induced effects can distort the low frequency part of the radio-source spectra even when tt is small (Sunyaev 
1971) and can influence the heating of electrons near active galactic nuclei and pulsars (Levich & Sunyaev 1971). 

The process of multiple scattering of radiation is described by a kinetic equation. This equation can be easily 
written if we neglect the induced scattering and/or the polarization (see, e.g. Nagirner & Poutanen 1994 and references 
therein). Induced scattering leads to the appearance of nonlinear terms in the kinetic equation. Since the polarization 
of radiation is described by four parameters, one must write a vector kinetic equation, i.e. a system of four equations. 
The kinetic equation combining these two effects (polarization and induced scattering) was not derived self-consistently 
up to now. The aim of the present paper is to fill up this gap. 

The first kinetic equation for Compton scattering with induced scattering was written by Kompaneets (1956) and is 
known under his name. He considered multiple scattering of homogeneous and isotropic radiation in infinite space filled 
with homogeneous, nondegenerate thermal electron gas. The gas was assumed to be non-relativistic {k^Tc <C mc^), 
the radiation to be rather soft {hv <C mc^), and the intensity of radiation to be a sufficiently smooth function of 
frequency. Because of small changes of photon frequency in a single scattering, the scattering integral was transformed 
to a differential operator by the Fokker — Planck method. The equation was rediscovered by Weymann (1966). More 
general Fokker — Planck equations were deduced by Cooper (1971) and Barbosa (1982) for more energetic electrons 
and photons, and by Molodtsov (1994) for the anisotropic, degenerate and moving electrons. For very cold electron 
gas, the Fokker — Planck equation was derived by Ross et al. (1978). 
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The induced scattering effects strongly depend on the angular distribution of the radiation field. Babuel-Peyrissac & 
Rouvillois (1969) generalized the Kompaneets equation to nonhomogeneous and anisotropic radiation. Their equation 
contains an integral operator in angular variables, but a differential operator in frequency. A simplified form of that 
equation was used by a number of authors to study induced effects in pulsar magnetospheres and radio sources in 
active galactic nuclei (Wilson & Rees 1978; Coppi, Blandford, & Rees 1993; Sincell & Krolik 1994; Lyubarskii & 
Petrova 1996; Sincell & Coppi 1996). 

As it is mentioned above, the induced scattering makes equations nonlinear. This effect can be simply taken into 
account in the equation if we limit ourselves to one parameter, intensity /, ignoring polarization. The scattering 
rate is increased by stimulated scattering by a factor 1 + h'^I{k)/2ck^ (where k is the photon momentum after 
scattering) comparing with the rate when only spontaneous scattering is considered. It is this factor that the authors 
of the mentioned works (Kompaneets 1956; Babuel-Peyrissac & Rouvillois 1969; Molodtsov 1994) have introduced 
into the equations. Often kinetic equations (e.g. radiative transfer equation) are deduced hcuristically from intuitive 
considerations and previous experience. However, in the case when polarization including induced effects is considered, 
our intuition and experience are not enough and the phenomenological approach fails. Therefore, one has to resort 
to deductive methods. Note that in nearly all the works on Compton scattering the electrons were assumed to be 
unpolarized and isotropic (or monodirectional). 

The first attempt to find out a rule for consideration of induced processes simultaneously with polarization was 
undertaken by Pomraning (1974). In terms of plane electromagnetic waves, he considered a source of polarized emission 
interacting with already existing radiation field and deduced a matrix describing the transformation of the Stokes 
parameters. The Stokes vector corresponding to the induced process was presented as a product of a 4 x 4 matrix 
N with the Stokes vector of the incoming radiation. Calculating the elements of the matrix N, he averaged the 
product of the quantities corresponding to the emitted (source photons) and passing by (field photons) waves over 
the distributions of their phases assuming that no correlation exists between the source and the field. This result is 
not applicable to scattering because the scattered radiation docs not have random phases. In the resulting expression 
for N, only the elements in the first row and the first column (i.e. related to the intensity) are correct in the case of 
scattering (Stark 1981). 

Derivation of the kinetic equation for the induced scattering of polarized radiation by non-relativistic electrons was 
given by Wilson (1978). He used the Maxwell equations, expansion electromagnetic field on set of harmonic oscillators, 
second quantization method and the non-relativistic perturbation theory. The scattering (a process of the fourth order 
in the approximation used) was represented as consequent absorption and emission. The products of the probabilities 
of these processes were replaced by the elements of the scattering matrix deduced according to the rules of quantum 
electrodynamics. While such a deduction is not self-consistent, the resulting equation is correct. 

Wilson's results were used by Stark (1981) to deduce the generalization of the Babuel-Peyrissac — Rouvillois (1969) 
equation for stimulated scattering to linearly polarized radiation (two Stokes parameters). Hansen & Lilje (1999) have 
corrected inaccuracies and misprints in Stark's work. Using a version of Wilson's equation containing only non-linear 
terms, Wilson (1982) and Coppi et al. (1993) showed that the induced scattering can enhance the polarization (compar- 
ing to the spontaneous scattering) producing also its strong frequency dependence. Induced Compton backscattering, 
for example, can amplify an incident ~ 1 per cent polarization up to ~ 50 per cent (Sincell & Coppi 1996). Accounting 
for the polarization of radiation may lead to the results qualitatively different from those obtained neglecting polariza- 
tion. A small change in the scattering rate (affected by polarization) can be strongly amplified due to the non-linearity 
of the process. 

The relativistic kinetic equation for the Stokes parameters taking into account the induced scattering was presented 
by Nagirner (1994) without deduction. The work of Wilson was not accessible for the author at that time. Shortly 
after that, the work of loffe made at the beginning of the 1950-ies was published (loffe 1994). In this work the equation 
was dcducicid (as loffe pointed out, with an approximate and a simplified method) for the so called polarization tensor 
which differs from the density matrix because the unphysical scalar and longitudinal photons are not excluded. 

In this paper, we derive the relativistic kinetic equation for the photon gas interacting with the electron gas. The 
way of deduction is based on the methods of Bogoliubov (Bogoliubov & Gurov 1947) and Baranger (1958). The scheme 
closest to ours is used in Silin (1971), where the non-relativistic kinetic equation for interacting electrons was derived. 
In our paper, at every step of the derivation all the expressions and relations are demonstrated to be either explicitly 
relativistically covariant or they can be changed to relativistic forms. The final equations are explicitly relativistic. In 
order to elucidate these facts and to introduce suitable for our purposes notations, we are urged to present a number of 
known relations and schemes. We thus fully describe the method of the derivation of the equations in a self-contained 
way. 

We make assumptions usual for the kinetic theory. Only binary collisions are accounted for. We assume the 
molecular chaos, i.e. the gas states are characterized by one-particle distribution functions. The characteristic temporal 

and space scales of a single interaction are assumed to be much smaller than the scales of significant changing of 
radiation and electron gas states. We use the principle of weakening of the correlations, which says that the correlations 
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between particle states are expressed in terms of the same one-particle distribution functions. There are no other 
limitations on the states of interacting particles. Simultaneously with the kinetic equations, we obtain the expressions 
for the scattering amplitudes and (well known) cross-sections. 

We first assume that both electrons and photons are polarized, then we average the equations over electron 
polarizations. An arbitrary anisotropy and nonhomogeneity of the radiation field and electron distribution are allowed. 
We also derive the kinetic equation describing the evolution of the electron gas. 



2. Density Matrix of the Photon-Electron Gas 

2.1. Operators of Creation and Annihilation of Photons 

According to the method of second quantization (see e.g. Bogoliubov & Shirkov 1959; Schweber 1961), the two types 
of vector-operators arc introduced satisfying the following commutation relations 

[af,{k),a,y{k')] = a^{k)a^{k') - a^{k')af,{k) = -g^^^kSik ~ k'), fi, i/ ^ 0, 1,2,3, (1) 

where k is the photon momentum, k — \k\ and the metrics {5^1^} = diag{l, —1, —1, —1}. The commutators of all other 
components are equal to zero, i.e. [a^(fe), a,y(A;')] = \aa{k),a^{k')] = 0. The commutator of the zeroth components of 
the vector-operator, ao{k) and ao(fc'), in equation (|l|) differs from others. Operators a and a are interpreted as the 
operators of annihilation and creation of a photon, respectively. Their product a(^)(fe)a(;^')(fe) (no summation) is the 
operator of the number of photons of momentum k and polarization A. 

Generally speaking, one has to quantize the photon field in a finite volume V (box), for example, of the paral- 
lelepiped shape with the sides equal to L^, Ly, L^, V = Lx Ly L^. Then in every dimension of this volume there can 
be a countable number of standing waves of the form (27r?i)~^/^ exp(ifcr/?i), where k — 2TTh{n^/L^,ny/Ly,nz/Lz) 
and nx,?T.y,nz are integers. Correspondently, the number of creation and annihilation operators should be countable 
and the commutation rules should contain discrete ^-functions, i.e. the Kronecker symbols. However, in the limit 
Lx, ijy,-£'z ~^ 00 the countable set of standing waves transforms to the continuum and all relations for the finite 
volume and its limiting case can be written in a unified way (Bogoliubov & Shirkov 1959). 

For example, the orthogonality condition for the discrete case with the integration over the box has the same form 
as in the continuum case, if the (5-function and the volume element of the three-dimensional momentum space are 
taken in the discrete form: 

J exp{z(fe - k')r/h}d^r - 5(fc - k'), 6{k - k') = ^-1^5„^„, 5„^„, J„^„, , d^fc = (2) 

Therefore, we keep nearly all the quantities and relations in the form corresponding to the continuum, while they are 
valid for the discrete case as well (i.e., for a finite box). 

Let us expand vectors a{k) — {ao(A;), a(fc)} = {a^^ik)} and gik) along the unit vectors of the four-dimensional 
basis whose vectors are orthonormal and their system is full: 

a{k) = a^^\k)e(^s^y a(fc) = a^^' (A;)e(;,) , e^x)e{x') = e(A)^e^^,) = gxy , e(A)^e[,^) = (3) 

(hereafter, summation over repeated indices is assumed). Raising and lowering the indices in brackets numbering 
the unit vectors are made in the same way as for the case of ordinary indices numbering the vector components. 
The coefficients of the expansion, a^'^\k) and a'^'^^k), are Lorentz invariants (i.e., they are scalars). They satisfy the 
commutation relations which follow from equations (^ and (^) 

[a^^\k),a^^'\k')]=-g^^'k5{k-k'), A, A' = 0, 1, 2, 3. (4) 

In order for the selected basis to become a polarization basis, their unit vectors should satisfy the additional relations 
that the two unit vectors corresponding to the transverse polarization should be orthogonal to the photon momentum 

fce(i)(fc) = fce(2)(fc) (5) 

where k = {k,k} is the photon four-momentum and the scalar product ke — k^e^ = keg — fc • e. As a result of 
such orthogonality, the photon momentum can be represented as a linear combination of the two unit vectors of the 
polarization basis 

fc=[Ae(o)(fc)][e(o)(fc) + e(3)(fc)]. (6) 
Then the following equality holds 

fca(fe) = [fce(o)(fc)][a("Hfc) - a^'Hk)], (7) 
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and similarly for the conjugate operator. Obviously, the unit vectors of the polarization basis are not arbitrary but 
related to the vector of the photon momentum k. However, the commutation relations are still valid. 

In classical electrodynamics, the transversality relation was given by equality to zero of the quantity (not an 
operator) formally coinciding with (|^). The equality of the amplitudes of scalar and longitudinal potentials (or equality 
of both to zero) was the consequence of that relation. In quantum electrodynamics, this relation is too stiff. One cannot 
require equality of a*^"^ (fc) and a'^-* (fe), since they satisfy different commutation relations which easily can be seen from 
relation (^. The solution of that problem was given by Gupta (1950) and Bleuler (1950) (see also Bogoliubov & 
Shirkov 1959; Schweber 1961). We use their method. 

Only the operators with indices 1 and 2 are physical and only those should enter the quantities having physical 
meaning. These operators satisfy the commutation relations 

[a^'\k)M''\k')] ^ [a(,)(fc),a(,,)(fe')] =<5,,.M(fc-fe'), s,s' = 1,2. (8) 



2.2. Photon States 

The states of the electromagnetic field in the second-quantized representation are described by the vectors with the 
arguments equal to given momenta and polarization projections of the definite number of photons. The initial vector 
that is used to obtain all other vectors is the state vector of the photon vacuum 'Sq. It is normalized and satisfies the 
condition of absence of photons: ^'o^o = 1, a^^\k)"^o = 0. The second equality is valid for any momentum k and any 
polarization A. 

All other photon state vectors can be obtained from the vacuum state vector by applying the creation operators 
but the vectors containing information about the non-physical scalar and longitudinal photons should be excluded. 
In the physical states, the operators of photon creation of these photons should appear only in the combination 
a*-°-'(fc) +a^^\k), so that there always exist equal amount of scalar and longitudinal photons. With the action of 
the operator a'-''-'(fc) — a'^^'>{k) on to such states, such photons mutually cancel out. The product of the operators of 
creation and annihilation in such a combination gives a^°^ (fc)a*^°^ (fc) — a''^\k)a^^\k) , which acting on to the vector of 
the physical states gives zero. The full product of the operators of creation and annihilation is reduced to a(fc)a(fc) = 
— 0(4,) (A;)a(s) (fc). Here the summation is done over s = 1, 2 (in the several following sections we use only lower indices 
s). We omit unphysical photons hereafter from our considerations. 

The physically allowed states of N photons with fixed momenta and polarizations can be represented in the form 

*si...s„(fci, fcjv) = ■-^a(s,)(fci)...a(s„)(fcAr)1'o- (9) 

The related arguments s and k can be moved simultaneously and rearranged in any order, since operators commutate. 

The orthogonality condition for the vectors of physical states follows from the commutation relations (||) and takes 
the form 

v[/,;...,^(fc'i,...,fc^)vI/.,...,„(fci,...,M = ^ E Il^k.S.^^'.Mkj - (10) 

ai ,...,Q!iv j=l 

where the sum is taken over all permutations ai, ...aN of indices 1, 2, N. These indices can alternatively be assigned 
to the arguments without primes, which is equivalent to the change from the sum over rows to the sum over columns. 
The results of the action of the creation and annihilation operators on to the state vector are given by 

1 ^ 

a(s){k)'ifsi...SM{kl,...,kN) = ~7^^ [fc(5ss,(5(fc-fc,.)] *si...s,_is,+ i...s„(fcl, ^r-l, ^r+l, few), (11) 

^ r—1 

a(^)(fc)«'si...s„(fci,...,/cjv) = %/7VTT 1'ssi...s„(fc,fei,...,fcjv). (12) 
The product of these operators acts on to the vector as follows 

N 

a(s)(fc)a(s')(fc')*si...SN(fcl' ■■■^kN)=^k'Ss'sr^{k'-kr)'^si...s^^iss^+i...SNif^lj ■■■,kr-i,k,kr+l, ...,fcjv). (13) 

r=l 

Arbitrary physical state with N transverse photons can be described by the vector 

—, C3i...s„(fcl, ...,fcAr)*si...s„('Sl: ••■,fcjv), (14) 

ki kisr 

where the sum is taken over the repeated indices Sj — 1,2. The coefficients of the expansion can be found as usual 
Fourier coefhcients. We assume that vectors of the physical state ( |l^ ) are dimensionless. They are normalized as 
^'w^jv = 1- In order to prove that, the orthogonality condition ( p^ should be used and the integral over momenta 
k[ must be taken applying the permutation property for indicies and momenta. 
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2.3. Quantization of the Electron-Positron Field 

In the second-quantization of the electron and positron fields (as in the case of the electromagnetic field) , the operators 
of annihilation and creation satisfy the relations of the form 

{br{p),blip')} = hr{p)hl,{p')+bl,{p')hr{p) = Srr'PoS{p - p'), for electrons, (15) 

and the same for the positron operators d^{p). Relation ( [l5| ) (and a similar one for positrons) contains anticommutators 
(not commutators), which is a consequence of the fact that these particles are fermions and satisfy the Pauli exclusion 
principle. All other binary combinations of the operators anticommutate. 

The electron-positron vacuum is defined with normalized vector $o as a state where there are no particles, i.e. 

6,(p)$o = 0, d,(p)$o = 0, $S$o-l. (16) 

Then the state with N- electrons can be defined via vector 

*ri...r„_(Pl,...,PArJ = -^=h\^{p^) . . .h\^ _{p (17) 

The order of the indices and corresponding momenta can be changed arbitrarily, but since the operators now anti- 
commutate, the sign of the vector changes depending whether a permutation is odd or even. The positron states can 
be defined in a similar manner. Since the scattering of photons by positrons does not substantially differ from the 
electron case, we can assume that it is only electrons that interact with photons and we do not introduce positron 
states. 

The operators of creation and annihilation act on to the state vectors in the following way 

1 

^r(p)«'ri...r«_(Pl,.-.,PAr_)=^^^(-l)^~^^rr,P0<5(p-Pj)$ri...T,_ir, + i...r„_(Pl,-.-,Pj_lPj+ (18) 
^r(p)^'Ti...r„_(Pl,-.,PAr_) = ^/nTTI ^ rr^. . .tm _{Pi Pli ■ ■ ■ ^P N J ■ (19) 

The products of operators give the following 

&r(p)&T'(p')*Ti...r„_(Pl,--,PAr_) ^^5r'T,Po5{p' - Pj)'^r^...T,-iTT, + ^...Tr,_{Pl^ -^Po-l^P^Pj + l^ -^PN-)- (20) 

Normalization of the state vectors is defined in a manner similar to photons. One has to change the order of the 
creation and annihilation operators many times in order for the former to be in the beginning and for the later to be 
in the end. The only difference is that a minus sign appears for odd permutation. The result can be written in the 
determinant form. For electrons, 

$|:,...^^_(p'l,...,p'ArJ*ri...r„_(Pl,...,PjV_) ^ ^\ (^r.rjP.0'5(p, - Pj )) • (21) 

The vector of an arbitrary state with iV_ electrons can be represented via an expansion over the elementary state 
vectors 

*JV- / ••■ Cri...r„_(Pi,...,PAr_)«'ri...r„_(Pi,...,PAr_). (22) 

J Pio PN-0 

The interchange of the arguments of the coefficients c changes their sign. The normalization of the state vectors, 
4>]y_$jv_ — 1, can be easily obtained from equations ( ^l|) and (^). 

2.4. Density Matrices for Photons and Electrons 

The density matrix of A^-photon state {N is fixed) can be defined as an averaged dyad product of the state vector 
with its conjugate 



^31,/ ^3t,. ^3i,,, _ 

(<;... s'„(fci>--->fcw)c.i....«(fei,...,fejv))«'.i....„(fci,...,few)«'.;....'„(fc'i,..., few)- (23) 



This product is an operator, since the elementary state vectors contain the creation and annihilation operators 
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Let us introduce the notation for the kernel of the density matrix operator 



(c:,...,^(fc;,...,fc'^K,....„(fci,...,fcA,)). (24) 

This kernel, being the A^-particle photon distribution function in the momentum representation, cannot be represented 
as a product after averaging of the expansion coefficients c. The density matrix is then written in the form 

and the kernel can be represented through the operator 

Only diagonal elements (i.e., elements with the same primed and non-primed indices and arguments) are real and 
non-negative. 

It is obvious from definition (^), that the density matrix is a self-conjugate operator, and it follows from (|2^) and 
( p5| ) that its kernel is a self-conjugate matrix, i.e. its Hermitian conjugation (complex conjugation and the replacement 
of the lower and the upper arguments with each other) is equal to the same kernel. The trace of the density matrix 
SppAT = (^Tv'I'Ar) = 1 by virtue of the normalization of the state vector. 

The density matrix carries a lot of information, much more than the distribution function or even the polarization 
matrix which can be expressed through the density matrix. However, in order to derive the kinetic equation, the 
distribution functions are not enough and one has to introduce a set of functions containing the groups of variables 
with the dimension less than N , but larger than 1. We call them truncated density matrices. 

Together with the iV-particle density matrix we introduce matrices whose kernels can be expressed through the 
integrals of the kernel of the original matrix. A kernel of order I 

s\...s[ ( k[ ... k'i\ _ 1 f d^ki+i d^kN f k[ ... k'l ki+i ... fcjv\ 

P'^-" Ui - kj- {N- ly. J - k^ [k, ... ki ... kr,) ^'^^ 

corresponds to the truncated Z-particle matrix. Integrals ( p7| ) are normalized according to normalization of the matrix 
Pat, so that 

Specifically, one-particle photon density matrix 

k'\ 1 /■ d^fcs d^kN fk' k2...kN\ ^^g) 



kJ {N-iy.J k2 kN \k k2 ...kN 

is normalized to the number of particles 
d^fc ,/fe 



/ 



The action of the product of the annihilation and creation operators on to the density matrix is reduced to 

N 



X '^si...Sj..isSj + i...SNi^lT--:^j-ljf'j^j + l^ ■■■i^N)'^s\...s'^{k'i, ...,k'j^)psl...s%[ 7^^)- (^1) 



1 ... KN 



Therefore, we get 

^sl^snC^Ii •■■'fcjv)a(s)(fc)a(s')('^') P '^s[...s'^{k'i, ...,k'^) 



1 ^ ' ' ' ' I I'h' h' k' h' 



j=l ^ J J ' 



H.,^(fc-fc,)p::::::-:---;;;-f ■•• - J . (32) 



N 
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The density matrices for particles (electrons) are introduced in a similar manner, therefore we immediately write 
down the expression for the joint density matrix for electrons and photons. It can be expressed through the kernel in 
terms of the following integral 



1 



d^fcid^fci d^fcAf d=^fc^ d^pi d^p[ d^PN^ d-^p'jv 



Pio Pia 

I I^n)^ti...t 



-Psi. 



PN-0 Pn^O 

iPi, --.Pn J'^l. 



.Sm ■Ti...Tm_ 



fel 



few 



{p'i,--,p'nJ- 



P'l 
Pi 



Pn^ 



(33) 



Changing the order of the photon arguments does not change anything, while changing the order of electron 
arguments in the state vectors or in the kernel of the density matrix changes the sign. However, if one changes 
simultaneously arguments in the state vectors as well as in the kernel, nothing changes since this procedure corresponds 
to just a change of notations. The sign does not change if one interchange lower and upper arguments, since the sign 
here changes even number of times. The truncated density matrix of electrons and their kernels can be introduced 
following the same procedure as in the case of photons. All the relations for the photon matrices are valid for the 
electrons too. 



2.5. One-particle Distribution Functions 

The kinetic equation for photons which we wish to deduce should be formulated for the one-particle polarization 
matrix depending on spatial coordinates, time, and photon momentum. Therefore, one has to make a transformation 
to such a matrix. We briefly describe the scheme for this transformation on the example of spinless particles. 

Up to now we used the momentum representation where the density matrices depend on a double set of momenta 
and polarization indices. The spinless states could be described in the coordinate representation. Instead of the matrix 
with two momenta arguments, we would have a matrix with two sets of space coordinates. The transition from one 
representation to another can be performed by Fourier transform. For the one-particle functions such a transition is 
defined by the formula 



(34) 



(35) 



where r = {ci, r}. The inverse transform is 

The two factors under the root are introduced since it is convenient when transforming to our notations. In the kinetic 
theory, the transition from the matrix with two space arguments to the usual distribution function is done via the 
Wigner function (Wigner 1932; see also Silin 1971), which is defined as follows 

p{p, r)= J d\e'P'"/^p,{r + v/2, r - v/2). (36) 

The Wigner function can be also expressed through the matrix with two momentum arguments 
1 f d^pidVl )r/h Pi+P'i\ J p'l 



The inverse transform reads 
'P'A _ VpoiP'oi 



PiJ {2T:hf 



J d'pd^re^^P^-P>/H (^p-Pl±Pl^ p{p, 



r). (38) 



According to the assumption of the small scale of the interaction between photons and electrons as compared with 
the macroscopic scale where the distribution function in equation (^8|) changes significantly, one can assume that the 
Wigner function does not depend on the space coordinates. Then both integrals in equation ( |38| ) can be taken and 
the matrix becomes diagonal in momenta 



p[l[)=PoM^Pi)p{Pi)- (39) 

Here we deduced these expression in three-dimensional coordinate and momentum space, i.e. non-relativistically. In 
order to introduce the relativistic generalization of the Wigner function (see, e.g., de Groot, van Leeuwen, & van Weert 
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1980), one has to define it for the generahzed momentum p which does not satisfy the relation = Pa ~ = rn^c^ ■ 
Then one can prove that outside of the surface defined by this relation in the four-dimensional momentum space, the 
function strongly oscillates (Zitterbewegung) and, on average, vanishes. Further consideration leads to equation (^9|). 

In the case of photons, situation becomes even more complicated, since the physical states with the transverse 
polarizations are connected to the momentum and, for the photon matrices with the fixed spin projections, transfor- 
mations of the type and (js^) are not possible. In that case, one has to return to the general polarization states 
which include non-physical photons. The elimination of the later should be done after the transition to equation similar 
to (|39|). Here we also meet with the difficulties to satisfy the relation fc^ = 0. 

Similar difficulty appears also in the case of electrons if one takes the solutions with a given spirality (i.e., a projec- 
tion of the spin on to the momentum) as the basic solutions of the Dirac equations. In that case, the transformations 
of the type (jsj) and ( ^7|) also are not possible. However, one can get solutions with the projections of the spin on to an 
arbitrary direction instead of the momentum. This can be achieved by the method proposed by Foldy & Wouthuysen 
(1950) (see also Schweber 1961). The mentioned transformations are then possible. The transition to the equation of 
type (|3^) is done in the same way as for the spinless particles and the spin indices are simply added. Without making 
these derivation, we take that 

"^'(pi) =poi'^(fi-pi)'°^'(pi)' (40) 

where the indices characterize the spirality. For photons wc make a similar approximation 

=fcj(fe'-fc)p:'(fe), (41) 

The matrices (p) and (fc) depend, of course, on time and space coordinates, but this macroscopic dependence is 
unimportant for the description of a single scattering act. 

In the normalization condition ( |30| ) the momenta k and k' are equal to each other. Formally, for equal momenta 
the (5-function in equation (|4l| ) becomes infinite. However, since we operate with the quantities in the finite volume 
(basic box), according to one has to take 

Then the normalization condition ( ^ ) reads 

/?4')^^^(^/'''^^^'^' (w/^^^^^(^)-y- 

When simultaneously V oo and N ^ oo while their ratio is constant, relation ( ^ ) transforms to the usual 
condition normalizing the polarization matrix. The sum p^{k) = p\{k) + P2ik) is equal to the double mean occupation 
number of the photon states 2p{k). This agrees with the fact that unpolarized radiation is described by the matrix 

p:'(fc) = (5.,.p(fc), (44) 
so that pl{k) = 2p{k). 

The normalization condition for the electron polarization matrix can be obtained from the equations analogous to 
(|^) and (|). The later takes the form S{p - p) = V/{2T:hf, so that 

^ U'pplip)^^. (45) 



The unpolarized electrons are characterized by the diagonal matrix 

P^r'iP)-^^^Sl'Mp), (46) 
where the distribution function in the comoving frame is normalized to the electron number density 



/ 



d'pfeip) = (47) 



We use the same notation, p, for both photon and electron matrices, but they differ by their arguments and indices. 
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In spite of the fact that equations ( ^ ) and (^) are deduced in a non-relativistic way they are relativistically 
covariant. Normahzation conditions (^) and ( ^5[ ) are written in the reference frames where the photon and electron 
gases are at rest. The corresponding relativistically covariant equations in an arbitrary frames are 



(2 



(48) 



where (fc) and (p) are the average photon and electron momenta in a given frame. 

Let us note that the interaction time Tq is also related to the (5-function. This relation is the relativistic counterpart 
of (0), i.e. its time- like form 



5{k^k) 



(49) 



3. Interaction between Photons and Electrons 

3.1. Equation for the Density Matrix 

Density matrix ( ^3|) satisfies the quantum Liouville equation (Landau & Lifshitz 1977; Silin 1971; de Groot et al. 1980) 

(50) 



= Hit)p{t) pmit) 



where H{t) is a Hamiltonian. This equation can be transformed to explicitly covariant form if we proceed from the 
Tomonaga-Schwinger equation (Bogoliubov & Shirkov 1959). In perturbation theory, the solution of this equation can 
be written as a series. It is shown in the text books on quantum electrodynamics that a part of this series corresponding 
to the scattering of a photon by an electron can be represented in the following form (index CE stands for Compton 
scattering by an Electron) 



C/^E(i,io) = Jdt' J dt" J d\' j dV"5CE(^,^^.) ^ J^y dV"5Ci=(r',r"), 

to to V V V V 

where V is part of Minkowsky space cto< < ct, ~Lx/2< x< £z/2, — Ly/2< y< Ly/2,—Lz/2< z< Lz/2, and 

Po p'o k k' 



cCE/ / //\ 



{2T:hf 
In the last expression 



-6t,(p')?^r(i>)a(,,)(fc')a(,)(fc)AA/;' 





p' 


-) 




p 


r" J 



(51) 



(52) 



^ S T 



and 



M 



S T 
S T 



M 



S T 
S T 



(53) 



dV^"'(p') 



ik') 



mc + Pv 



—V 



iO 



3(s)(fc)- 



mc + Pv 



iO 



3(,,)(fc')e'^--""--^/'' 



Jip'r'-pr")/h -ip (r'-r")/h 



^{Pl 



(54) 



Here v7{p), the elementary solutions of the Dirac equation describing electron, are the columns containing four 
functions, u'^{p) — [u'^{p)]^j'^ and p = pj = Po7° ^ P ■ 7 — Pfj.l'^, where 7^* are the Dirac matrices (e.g. Schweber 
1961; Berestetskii, Lifshitz, & Pitaevskii 1982). The scalar product kr — kct — k-r. The imaginary term iO determines 
the definite rule of circuit of singularity in the denominator. It is obvious that matrices ( |5^ ) and (^) are symmetric 
relative to r' and r". 

U V ^ 00, to ^ —00, t ^ 00, the electron Compton operator J7^^(— 00,00) = contains two integrals over 
four-dimensional space (r' and r"). Two four-dimensional ^-functions appear after calculating these integrals. One 
of them disappears when we take the integral over four-dimensional momentum of virtual electron p^, whereas the 
second one reflects the conservation laws. Thus integrals over the whole time-space (xq — ct) 



dV / d^r"Af: I 



,fk' 



dV /dV"Ai:';f 1^' 



: {2Trhf5{p' + k^-p-k)M'^'^ 



(55) 
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where the amphtude of the process 



= mcu^ (p') 



(fc')- 



(k) 



(k)- 



k' 



(fc') 



The amphtude has an obvious property 

k 



Then for matrix S'^^ the following expression is obtained (Bogoliubov & Shirkov 1959) 



^CE 



9 

e I 

he 2'K 



k k' 



Po Pa 



(56) 



(57) 



(58) 



It is easy to show that the operators Uy^{t,to) and 5*^^ being conjugated only change the sign. 

If we omit all terms which are higher than the second order and those which do not conserve the number and/or 
the quality of particles then the equation for the photon-electron density matrix describing Compton scattering by 
electrons takes the form 



dpjt'] 
dt' 



At' = p{t) - p(io) - C/^^(i, to)p(io) + pihW^""' 



it,to). 



(59) 



The derivative in equation ( p9[ ) describes all variations of the density matrix which have very different scales 
(hierarchy of scales). The duration and the characteristic scale of the interaction are much less than the scales of 
macroscopic changes of the matrix and this fact is the basis of the Bogoliubov method of deducing the kinetic 
equations. Thus we assume that the density matrix p changes negligibly during the interaction. Following Baranger 
(1958), we take the value of the derivative at time t out of the integral in the Ihs of equation (p9|), i.e. the Ihs is now 
Todp/dt. We write the full derivative because it describes now macroscopic changes. Thus we get the initial basic 
equation 



To 



dpjt) 
dt 



= p{t)-p{to) = U^''{t,to)pito) - p{tom''{t,to) 



(60) 



3.2. Equation for the Kernel of the Density Matrix 



Let us write equation l^v^) not in the operator form, but via the kernel of the photon-electron matrix keeping temporarily 
the number of components fixed (specifically, N photons and N- electrons). Simultaneously we substitute the time 
argument t of the density matrix instead of t^, because this argument is macroscopic time. 

According to equation ( p6| ) for photons and analogous for electrons, we multiply equation ( |60| ) by the corresponding 
state vectors and their conjugates on the left and right. We get 



ki ... kpf 
ki ... k^ 



Pi 



Pi -Pn^ 



,(fci,...,A;Ar)$t^ ,.^_(Pi,...,Pjv_) 



X [U^''{t,to)p{t) - pm^''{t,to)]^s[...s'Jk[, ...,k'j,)^r,,,,r,^_{p[, ...,p'^J^^^ 

g2 _ r 

:^s,...sr,iki...kN)<i>l,„r, (Pi-PnJ / dV' 



V V 

dXd^ d^fc^ d^fc^^ dV/ d^p'C d^PN_ d^PN 

h" 



d r 



k k' Po p'q 



/ ^ S T 





p' 


-] 




p 


r" J 



Ps 



■■^N ■ ' 1 

■ s'L,r'' 



N 



k" k" 



Pi" 
Pi' 



■Pn^ 
■P^_ 



yi yii yi yii ^ii ^iii ^ii ^iii 
«-! 'I'l '•'AT i^N PqI Poi PoN- PoN- 

hl,{p')br{p)a(s'){k')a^s){k)-^s'^...s^^^^^^ fc^)$,,,..,._ (p'/, ...,p^J^,,„...,.,(fci", fc^)*^,,,, (pi" 
(fci', ...,fc^)$......._ (pi', ...,p^jW,;....,.,(fei", ....k%)'^l,,, ^,,,_{p'i' , ...,p-J 



X b\.{p)hr'{p')a(s){k)a(s'){k') 



,{k'i, ...,fciv)*r(...T;,_ (pi, ■••,Piv_)- 



,PAr_ 



(61) 
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Using the normalization conditions ( [LOD and (21) for the state vectors as weh as formulae ( |13| ) and (gO|) describing the 
action of the creation and annihilation operators on to the state vectors, and integrating over all photon and electron 
momenta with two and three primes, we get 



To 



d s[...s'^,t[. 

^Psi...Sjv,Ti. 



N 

kN 



p'l 
Pi 



■ Pn^ 



N 

EE 

T-=l j = l 



kSss'S{k - k'^)p()SrT'S{p - p'.)ps\ 

" 3 



Sr+1- 



dv" 



d^fcd^fc' A^p dV 



Po 



p'o 



■'^ S T 





p' 


-) 




p 


r" J 



.Sk ,Ti...Tj-iT Tj-i-i ...rjv_ 



fcl 



k' 



k 



.. k 

p'l- 

Pi ■■ 



p' 

■Pj 



p'l ■ 
Pi ■ 

■P'n^ 
■Pn^ 



P 



■Pn^ 
■Pn^ 



(62) 



Deriving this equation, we used a number of times the possibility to change the summation order between rows and 
columns (see derivation of eq. jlOj) and as a result the factorials canceled out. The difference of two values of density 
matrix in equation (nO) gives the difference of corresponding matrix elements which we do not write. 



3.3. Equations for the One-photon Matrix and the Correlation Matrix 



Let us write now the equation for the kernel of the truncated photon matrix of the first order which follows from 
equation (|6^). We equate all corresponding upper and lower polarization indices and momenta of electrons and all 
(except the first) indices and momenta of photons, i.e. we sum over indices and integrate over momenta. The result, 
according to the definition (p9), is then divided by {N — 1)1 



^° AiP'^ Ui 



t \ 



{N-l)\N-\ 



N 



dV 



dV" 



d^fcd^fc' d^p dV ,^Yfc' 

~J S T 



k k' 



Po Po 



k'Ss'sJik' - fcl) XlP0'^^'rj'^(p' -Pj)/0s's2.'.'.'s«.Ti\' 



i=i 



^k'Ss's,.Sik'~kr)^PoSr'TjS{p'-pj)pll 



S2- 



kS,ss{S{k - k[)^poSrTjS{p - p'j)pl^T2.''X] 



.sjv,ri. 



■ rN_ 

.Tjv_ 



.rjv_ 
.rjv_ 



fc'fc2 
fclfc2 



■ Tn_ 
,.rjv_ 



fc'lfc2 
fe fc2 



k[k2 
k fc2 

fc 



d^A;2 d^kN d^Pi d^p_ 



kN POI " PON- 



kN 
kN 



Pi 
Pi 



Pi 

p 



■Pn^ 
■Pn^ 



kN 
kN 



Pi 
Pi 



Pj 

p 



Pn^ 
Pn. 



N 



^kSss,.S{k - kr)y ^Po5tt,5{p - p^)ps\ 

r=2 j=l 



SlS2- 
S2- 



■ Sm,ti...t ...tn_ 

■ Sj^.Ti...Tj...TN_ 



kN 

kN 
k[k2 

fclfc2 



Pi 
Pi 



P' 
Pj 



Pn^ 
Pn^ 



fc' 



kN 
kN 



Pi 
Pi 



P' 
Pj 



Pn^ 
Pn^ 



(63) 



Here in summation over r we separately write the first term corresponding to r = 1. Since every integral in the sums 
over r and over j gives the same result, the summation is reduced to the multiplication by the number of terms in 
the sums, so that the factorials cancel out. The sums over r = 2, N give the same kernel of the truncated matrix 
Psjs^'r and cancel out. The resulting equation in the right hand side contains only kernel of the matrix depending on 
variables of one photon and one electron 



^°di^- Ui 



dV / dv 



d^fcd^fc' d^p dV 



k'Ss'siS{k' ~ ki)pl]'l 



k k' 



Po Po 





p' 


-) 




p 


r" J 



ks,s[Sik-k[)pi-:-; 



(64) 



Let us now derive the equation for plUrl- We follow the same procedure as when deriving equation (|6j), fixing 
the characteristics (momenta and indices) of one photon and one electron in equation (^2|). One has to separate these 
variables when we consider the action of annihilation and creation operators on them. The resulting expression contains 
four couples of terms with pluses and minuses. The first couple does not have summation at all (i.e., corresponds to 
r = 1, j = 1), the second and the third couples contain summation over the photon and electron variables, respectively, 
and the fourth couple contains summations over variables of both interacting particles starting from r = 2 and j = 2. 
The terms in the last couple (with double sums) give the same function Ps\s\t\t and annihilate. We take as the Ihs 
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of the equation the difference which stands on the second place of equation (Q) . We omit here the argument t of the 
density matrices in the rhs of the equation based on the assumption formulated in the next section. We get 



P'l 
Pi 



tn \ = i 



dV 



dV" 



d=^fc d^fc' d?p d?p' 
k k' po p(, 





p' 


-) 




p 


l" J 



k'5,.s,5{k' - ki)p'oSr'rAp'-Pl)Ps':"M'''' 



I if k' 

k5,s[S{k - k[)poSrr'^S{p - p[)pl^i'!;.i^^ 



k'6sisAI^'~k,)pt^:tr'^''' 



p[p' 



Pi 



p'lp' 

Pi p 



fk[ 


k' 


")] 




k 


PiJ. 



(65) 



In order to proceed further we will use one additional assumption described in the following section. Then we explore 
once more the smallness of the interaction scales in comparison with the scales of macroscopic changes. 



3.4. Approximation of Weal<ening of Correlations 



The matrices depending on variables of a number of particles contain the information about correlations between these 
particles. In kinetic theory, usually one applies an approximation of molecular chaos (Silin 1971). According to that 
approximation, there are no correlations before the interaction, i.e. 



Psi.Tl 



p'l 
Pi 



to 



Ps] 



'.(k'l 
ki 



(66) 



An interaction creates the correlations. However, according to the principle of weakening of correlations which is 
satisfied for sufficiently rarefied gases, the correlations are accounted for only in equation (|6^ ) for the one-particle 
photon matrix via kernel p entering the right hand side of the aforementioned equation. This kernel is assumed to 
characterize the electron and photon states after the interaction. It can be represented via the same kernel before the 
interaction and the correlation function as it is done in equation (|65|). 

The two-particle truncated matrices can be presented as products of the one-particle density matrices accounting 
for the exchange effects (symmetrical and anti-symmetrical forms for photons and electrons, respectively) 



PsiS2 



ki k2 
ki k2 



Ps, 



-Ps\ 



Ps, 



''"l''"2 
PtiT2 



Pi P2 
Pi P2 



Pi 



Pt2\ 

\P2 



For matrices of higher order entering equation (pq) , similar equations hold 



Psisl.rl 



ki k2 
ki k2 



^ PS1S2 



ki k2 
ki k2 



Pri 



PS1,T1T2 



P'lP'2 
Pi P2 



= Ps 



PriT2 



p'lp'2 
Pi P2 



(68) 



It can be shown that these matrices satisfy all normalization conditions. 

Now we substitute equations (|66|)-(|68|) into the rhs of equation (65), and the resulting expression into equation 
(p4). At this step, after all substitutions, we will use the fact that the matrices in the rhs of equation (65) do not 



depend on time during the interaction. Then we can take the integrals over the space-time variables enter the matrices 
M only. Simultaneously we let V 00, —00, keeping the upper limit t finite. As a result of this procedure four 
integrals over space give four three-dimensional ^- functions as in (p5|), namely two pairs of the type 



6{p' + k' - p^)5{p + k-p^)^ S{p' + k' -p- k)5{p + k 



(69) 



and 



S{p' k - Pv)6{p k' - p^) = 5{p' + k' p k)5{p - k' - p^ 



(70) 



Both combinations give the (5-functions reflecting the momentum conservation. The second set of 5-functions gives a 
possibility to take the integrals over p^. 

Integrals over time with the upper limit t can be taken independently. They have the form 



J{k'+p'o-p°)ct'/n^^/ ^ ^i{k'+p'o-p°)ct/h 



irhSik' +p'„-p",)-- 



k'+p'o-p° 



(71) 
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Integrals over with the second terms of this type must be calculated as the Cauchy principal values. These terms 
do not ensure the energy conservation law. They vanish if the limit of integration t is equal to oo. 

The terms with absence of energy conservation in kinetic equations, as was noticed by Nagirner (1994), appear 
in Silin (1974) and Bomier (1991) who do not discuss this fact. In their papers, the interaction between electrons 
(Silin 1971) and between photons and atoms (Bomier 1991) does not depend on time. In our case, the Hamiltonian 
is more complicated and without the energy conservation law we cannot calculate the scattering amplitudes and the 
cross-sections of the process. 

We suppose that such terms are not physical and must be excluded. The interaction between photon and electron 
has finished before the density matrix changes noticeably and the integration limit t must be taken infinite. Then 
the resulting integral is given by equation (p5|). However, this procedure gives the cross-section twice as large as the 
correct value. To overcome this difficulty we continue the upper limit t to oo but take half of the result. 

The next step is to use for the last time the assumption of the small scale of the interaction. Instead of the one- 



particle kernels we use now their diagonal in momentum forms ( |39D and (41). Accounting for relation {]^, the Ihs of 
equation (p3) then takes the form 



^°dtP'^ Ui 



Toh6{k[ - ki)^Osl{ki) = ki6{k',-k,) — ^oi{k,). (72) 
at c at 



After all the substitutions mentioned above, the rhs of equation ( [6^) contains 22 terms. The two terms are the 
result of the substitution of (|6^). All other terms in equations (|^ and (|65|) come in couples, so that the total number 
of terms is multiple of 4. Every term in the first line in the square brackets of (^) gives two terms, and every term 
from other two pairs of terms gives four terms, resulting in 20 terms. 

The first two terms containing linearly the amplitude M given by equation ( p6| ) cancel out. This follows from the 
equality 



MV 



= Sss'SrT', (73) 



which can be easily proven using properties of the solutions of the Dirac equation. Equality ( [73| ) leads also to the 
cancellation of the eight terms arising from the first terms in the rhs of equations (|67|). 

All the remaining twelve terms contain eight three-dimensional integrals over photon and electron momenta. There 
are two four-dimensional and five three-dimensional 5-functions in the integrand. All three-dimensional ^-functions 
disappear when taking the integrals. The resulting expression contains three three-dimensional integrals (one over 
photon momenta and two over electron momenta), and two four-dimcnsional (5-functions, 6{p' + k' — p — k) and 
^ik!i ~ ki)- The first one reflects the energy and momentum conservation laws. The second one cancels out since it 
also appears in the Ihs of equation ( |6^ ) according to equation ([72|). 

In the final expression, four out of twelve terms contain products of one photon and one electron matrices corre- 
sponding to the spontaneous scattering. Two of these terms entering with the plus sign are responsible for the emission 
and the two terms entering with the minus sign are responsible for the attenuation of the photon beam due to scatter- 
ing. The remaining eight terms reflect the exchange effects. Four terms containing the product of one electron and two 
photon matrices correspond to the induced scattering, and four terms containing the product of one photon and two 
electron matrices reflect the exclusion principle. The induced scattering terms have the same signs as the spontaneous 
terms, while the exclusion principle terms have the opposite signs. After an elementary but very long deduction, we 
obtain the sought kinetic equation. 

4. Kinetic Equations for Photons and Electrons 

4.1. Kinetic Equation for Polarized Photons 

As we mentioned above, all the remaining terms in the rhs of equation (|6^ ) after substitution of equation (^6|) contain 
the factor (5(fc^ — ^i) which also exists in the Ihs of the equation. This factor therefore cancels out. We can now restore 
the dependence of the photon matrix on time and spatial coordinates. Let us also change the designations, writing the 
photon matrix in the Ihs in the form pl°{k). Rewriting the derivative over the line of sight as the full derivative, the 
Ihs of the equation takes the covariant form 

k^pi{k,r,t)=f~+k-v)pi{k,r,t). (74) 
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In the rhs of the equation we change the indices and add and subtract two pairs of identical terms containing products 
of two photon and two electron matrices for symmetry. The resulting kinetic equation takes the form 



kVpZ{k,r,t) 



2 {27:h)3 J k' po p'o 



d3fc'd3pd3y^^^,^^,_^_^^ 



5i,+pUk') 



st'+pi'ik') s;'-p;'ip') 



MX' 



^ ^"snr"' 



k' 

k 



pi„{k')pUp')Myr' 



s;„-pi„ip') 

k 



MV, 



Ap7(k)pl iP) 



k' 



St„+ptAk) 5l,„~pl,„{p) 



S T 

s" t' 



P 
P' 
k' 
k 



(75) 



Here vq = e"^ /mc^ is the classical electron radius. The arguments r and t of the matrices in the collision integral are 
omitted. The terms in the first two lines in braces describe the photons outgoing from the state with momentum fc, 
the other two lines describe the photons incoming into this state. In order to outline the matrix character of Kronecker 
symbol we use further the notation J| , where the lower index represents the rows and the upper index represents the 
columns. 

Note that the products of four matrices are added just for the symmetry of the equation. If these four terms are 
dropped and we change some indices, the equation takes the form 



kVpZ{k,r,t) 



,2^2 



2 (27r;i)3 



d3fc'd3pd3y^^^,^^,_^_^^ 



MX, 



k' Po P'o 

'Mi,;:, - Mf:;„p;'(p') + M!„^„fyi{k!) 

Ml„l„ - pl„{p')My,, + pi„{k')MI,;,?\^ p'hk)pf{p) 
2MtJr„ pi {k')pl {p')MX + 2Mf„;„ pi {k')pl {p')MX"p^r'" iP) 

M^„pi {k')pi (p')Mi'';"4 (fc) - pio miXr" pi Wr ip'WsK' ] 



^"so t'" 



(76) 



In the last equation we omit the arguments of scattering amplitudes M for brevity. These arguments are the same 
and in the same order as in equation ([zs]). 

We note, that a possibility to calculate the scattering amplitudes appears after we choose the polarization basis 
for photons and the axis where one projects the electron spin. As the axis one can use the electron momenta. The 
situation is more complex for the polarization bases. We will give detailed description of that later and now turn to 
the equation for the electron distribution function. 



4.2. Kinetic Equation for Polarized Electrons 



Now we also can write down the kinetic equation for electrons. It can be deduced following the same procedure as 
in the case of the photon equation. However, due to the symmetry of equation ([75|), we can write this equation by 
analogy. One just removes one integral over electron momentum and adds one integral over photon momentum as 
well as sums over photon polarization and removes summation over indices of electron polarization (with the external 
momentum, p). The resulting relativistic kinetic equation takes the form 



PYpro(p,r,t) = - 
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(77) 
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It is easy to check that in the arbitrary time moment, t, and in the arbitrary point in space, r, two relations are 
valid 

vjkpllik,r,t)^^0, vjpp-!;^ip,r,t)^ = 0. (78) 

These relations reflect the continuity equations for the photon and electron number and are the consequence of the 
conservation of the number and quality of particles during the scattering. It is also easy to check the validity of the 
energy-momentum conservation for Compton scattering 



= 0. (79) 



We note that conservation laws (|7g) and (|79| ) are valid only for the occupation numbers, but not for the polarizations. 
4.3. Scattering Amplitudes 



Matrix elements (^6|) for Compton scattering (amplitudes) can be rather easily calculated for a special choice of 
polarization bases — internal bases. We take the unit vectors of these bases in the following form. For the photon of 
momentum k 

ej" (fc) = ei"afc) ^ ^ {kk'p,pok xk' + kk' xp + k'pxk). (80) 

Here the dimensionless scalar products 

^ = kp/m^(? ^ It p' /inn?(? ^ S^' = It p / rn? (P' = k p' / rn^ (? ^ q = k It /rn^c^ , (81) 



and A = \/2^^'/q — 1. The double definitions of ^ and ^' in equation (^ij) as well as the relation S, = S,' + q are the 
consequence of the conservation laws. 

The values <^ and ^' are the photon energies before and after scattering in the frame where the electron before 
interaction is at rest. If we denote the cosine of scattering angle in this frame as /io then 



9 = ?C'(l-Mo), C' = T— tI -y A=J\±^, ^0 = 1+1 1 (82) 

l + C(l-Mo) V 1-Mo C C 

For the photon of momentum k' , the unit vectors are very similar 

efi) (fe') = 4i) (^) ' £(2) (fc') = -42) (fe) • (83) 

One can add a vector of the photon four- momentum k (or fc') multiplied by any real number to the polarization unit 
vectors. By doing so, one can make the time components of these unit vectors (as well as of the external vectors e(fe)) 
equal to zero. For such vectors, the scattering amplitudes (matrix elements) which are marked with a circle atop of 
the letter, are 





° 12 

= Mil = 




°21 


= ^mH 


. i 

^'-^ 


°ii 


°12 

= -Mil 


.c 



° 12 

Mil 


°ii 

— ^"12 














°22 

Mil 


°21 
= Ml2 = 




- U-W-), 


°21 

Mil 


°22 

= Mil = 




(84) 


°12 

Mil 


°11 

= Mil - 








°21 

= -Mil 







Here we denoted 



v± = Vr + iVr' + 1 ± Vr - iVr' - 1, w± = Vr + iVr' - 1 ± Vr- iVr'-i-i, 



4^A'(A' + Z) y g Vl-Mo 

2vC C/ r' 2^\£,' (,) r' ' ^ 

Note that the scattering amplitudes enter equations ( |75| ) and ( |77| ) without summation on their indices only with the 
products of four density matrices (two photon and two electron ones), but as we mentioned these products annihilate. 
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In all other terms, as we can sec in the equation for photons in the form ([76|), there are summations over at least one 
pair of indices of scattering amplitudes. 

If indeed the summation is over only one pair of indices then the resulting sum depends on the 6 remaining indices 
and contains 2^ = 64 terms altogether. If, on the other hand, the summation is fulfilled over two pairs of indices, then 
the sum contains 2^ = 16 terms. The term which describes the attenuation of radiation belongs to the last category 



''s't' — ^"s' r' 



(86) 



where the sum is taken over one pair of photon indices and over one pair of electron indices. Here one gets 



''11 — '■12 — ''21 ^ '■22 — 



, 11 _ :r 12 _ 

til — '■12 — 2 ^ 



° 21 _ ° 22 _ 
'■21 — ^22 — 2 ' 



'11 



' 12 



^21 



^22 



,21 
' 11 



12 
'22 



22 
'12 



^21 



(87) 



where B = ^/e' + eve 



4.4. Transformation Matrices 

We have deduced the formulae for the Compton scattering amplitudes (^) in a specially chosen internal polarization 
basis (pO|). However, the expressions ( ^^ and (56) for these amplitudes contain unit vectors of some external bases 
which should reflect the geometry of the scattering medium. The photon density matrices entering the kinetic equations 
([75|)-(|77|) are connected to the same (external) vectors. Therefore, we must find transformations between the external 
and the internal bases. 

Two sets of polarization vectors e(^g-){k) and e^'," (fc) together with the vector k/k form two three-dimensional bases 
with one common unit vector. Consequently, the external vectors and the internal vectors are connected to one another 
via relation 

e(.)(fc)-?,':(-x)4")(fe), 

where the angle x is defined in the following manner 

cosx = -e(i)(fc)e(']")(fc) = -e^2)ik)e(^)ik), sinx = e(i) (fe)e('2') (fc) 

The specific expressions for cosx ^'^d sinx are given in Nagirner 
defined for the plane-parallel medium. The transformation matrix 

cos X sin X 
— sin X cos X 

has the following properties 

^.':(x) = ^^.(x) = ^.^(-x) = r,(-x). (91) 

These properties are obvious if the numbers 1, 2 are substituted instead of indices. Matrices of this kind arc commutative 
and multiplicative, and their transposition is equivalent to the inversion or changing the sign of the argument. 
Multiplying both sides of equations (^8|) by the Dirac matrices we get 

e(,)(fc) = l-:.{~x)ei'^^{k), e(,o(fe') - (-xOe^) (^O- (92) 

The transformations of bases related to the photon momentum fc' are analogous to ( p8| ) and the angle x' is defined 
analogously to x in equation (|89|). Then, for amplitude (p6|), we have 



{ifix}} = 



-e(2)(fe)e(\")(fe). 
Poutanen (1993) where the external vectors are 

(90) 



S T 

s't' 



'5(5) (fe) 



mc + p - 



k' 



e(,/)(fc') - e(y)(fc') 



mc + p + I 



'■e{s}ik) 



Is''.{-X)ls'.i-X')^u^ ip') 



6(7) C^)- 



P ^'eP)(fc')-eP)(fc') 



mc 



p + k 
mcS, 



e(V")(fc) 



(93) 



Using equations ( |57| ) and (|9^) we obtain the following expression for the matrix element with the primed upper 
momenta 











p J 


p'j 





i''r'hx)i-:.i~x)Mrr 



^i''r'hxm-x) 



(94) 



These formulae must be substituted into equations (|75|) and ([7' 
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4.5. Equation for Photons Interacting with Unpolarized Electrons 
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If electrons arc unpolarized (i.e., equation is valid), then the matrix pl^ (p) must be replaced by the product 
(27r?i)^/e(p)/2. The factor — (p) is replaced by [1 — {2ttTi)^ fc{p)/2], and one can sum over the electron 
polarizations. For the external bases. 



-L „/ 



S S 

s' s' 



k' 



ui'„i(l, 



(95) 



Here, one sums over the repeated indices r, t' = 1,2. After summation over electron polarization the Compton 
scattering cross-sections of polarized radiation cannot be represented as a product anymore. This reflects the transition 
from the "pure" scattering of photons and electrons with fixed polarization states to the mixed one. 

The expressions for the cross-section s (p5[ ) have the simplest form in the internal polarization bases (^^ and (^3|) . 
The matrix elements are then given by (R4)and therefore 



5,s s" ( k 



o 



P 

P'J' 



(96) 



where the sums are taken over the repeated indices t, r' = 1,2. 

Only 8 out of 16 elements of matrix (^^ are not zeros, and only 5 elements are different 

Til ^ ^[B + 2 + 4(Mg - 1)], T^^ = -{B + 2), = = -{B - 2), 



nl2 
■ 12 



°91 1 °99 °11 1 

tH = -{B + 2)^.0, TiT = tH = -^{B- 2)/io. 



(97) 



The same expressions are obtained if they are derived directly for unpolarized electrons by means of traces and 
projection operators (see e.g. Berestetskii et al. 1982). 

Combination ( p3| ) and ( p^ ) with (^5|) gives the law of transformation of the cross-sections of polarized radiation 
by unpolarized electrons. One can write a transformation, for example, in the following way 



rps s"( ^ 



p' 



p' 



, )lrf'ixKs'-"(x')- 



(98) 



The transformation matrices can be written in various forms using equation (pi|). 

Keeping in mind the formulae expressing transformation laws and omitting the arguments of the T-matrices, we 
rewrite the kinetic equation for photons with unpolarized electrons in the form 



9 9 9 
Tn m C 



—0[p +k ~ p- k) 



k' Po P'o - 

{/e(p') [1 - i2nhfup)/2] 2T:jpiik') - Up) [1 - {2^hfup')/2\ [pi;,{k)r^,'^ + T:;;:„p^ik) 
+ [Up') Up)] [pi:;{k)Ty}pi {k') + T:/,,pi:„{k')pt"{k)] } . (99) 

This equation can be easily rewritten in a more customary for astrophysics form in terms of the Stokes parameters. 



4.6. Kinetic Equation in Terms of the Stokes Parameters 

The transformation from the polarization matrices to the Stokes parameters is done with the aid of the formulae 

where ni = pl{k)/2 = p{k) is the mean occupation number of photon states, uq, nu, uy are the corresponding 
characteristics of photon polarization which can be transformed to the standard Stokes parameters (having dimension 
of intensity) by multiplying them by the factor 2fc^c/(27r?i)^. The equation for polarized radiation interacting with 
unpolarized electrons takes the following form in terms of the Stokes parameters 

f d^k' d^n d^n' 

kYn{k,r,t) = -^m'c^ / — ^^5(p + fc - p' - fc') {/c(p')[l - (27r?i)Vc(p)/2]L(-x)FL(x')n(fc') 

+ [Up') - /c(p)]N(fe)L(-x)FL(x')n(fc') - /e(p)[l - {27Tn)^Up')/2][F + L{-x)AL{x)]n{k)} , (101) 
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where 



N(fc) = 



/ ni 
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ny \ 




ni 








nu 
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(102) 
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L(X) 



/I 







^ 

COS 2% sin2x 
- sin 2x cos 2% 
1/ 



(103) 



The matrices L(x) have the same properties as matrices (|90|). They are commutative and their transposition is 
equivalent to the inversion or changing the sign of the argument. Note that the same argument x appears with the 
opposite signs in matrices L (around matrix A) in the last term of equation ( 101 ). 

As we have mentioned above, the equation of the form (101) for non-degenerate electrons (not accounting for the 
exclusion principle) was given in Nagirner (1994). The term with matrix A was omitted there. It can be easily shown 
that this term disappears if the electrons have an isotropic distribution (even if they are degenerate). In the major 
part of Nagirner's paper, the isotropic electron distribution was assumed, so that this omission did not introduced 
any errors. The ma trix A was omitted also in papers of the authors (Nagirner & Poutanen 1993, 1994) where the 
rhs of equation (101) was averaged over the directions of electron momenta and the five functions describing the 
redistribution of radiation in frequency, angles and polarization states were obtained. However, these papers were 
devoted to the isotropic electrons as well. 

Electron isotropy means that the medi um i s locally isotropic, and therefore the attenuation (i.e. the terms corre- 
sponding to the last square brackets in eq. |101]) is described not by the matrix, but only by the scalar determined by 
the Klein — Nishina cross-section, F. The cross-section averaged over momentum of scattered photons can be found 
in Nagirner & Poutanen (1994). Hence we can conclude here that polarized electrons can infiuence both the linear 
and circular polarization by means of attenuation through scattering. Unpolarized, even non-isotropic, electrons can 
introduce and change in this process only linear polarization. 



5. Conclusions 



In this paper, we have deduced the relativistic kinetic equations which describe the behavior of the rarefied photon 
and electron (or positron) gases interacting with each other via Compton scattering. We accounted here for stimulated 
effects for the photons and for the exclusion principle for the electrons. We considered arbitrary polarization states 
of photons and electrons. We presented also the kinetic equation for polarized photons scattered by unpolarized 
electrons in terms of the Stokes parameters. The expressions for the scattering amplitudes and cross-sections are 
derived simultaneously. There are no limitations on photon and electron energies. 

Note that all the deductions were made by means of relativistic quantum electrodynamics methods and all the 
equations obtained are relativistically covariant. For particular scattering problems they, of course, must be adapted to 
the geometry and symmetries of a medium and initial and boundary conditions (see e.g. Nagirner & Poutanen 1993). 

Finally, we would like to notice that the factors 1 -I- p{k) and 1 — (27r?i)'^/e(p)/2 accounting for the induced effects 
for photons and the exclusion principle for electrons in the scalar kinetic equations, should be replaced by the new 
factors (5f + (k) and (5^ — (p) in the kinetic equations for polarized photons and electrons. This rule can help to 
formulate kinetic equations for more complex processes for which a direct derivation could be very complicated. The 
processes in strong magnetic fields can serve as such examples. 
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^ Let us not e he re that the generalization of the Babuel-Peyrissac — Rouvillois equation to polarized radiation (i.e. a limiting 
version of eq. [101] for weakly relativistic Maxwellian electrons and soft photons) was deduced in Nagirner (1994). The English 
version of the paper contains many misprints in the generalized equation. In the Russian version, n should be substituted instead 
of n in the first row of his equation (41) and the brackets b efore n" and after h'l in the third row should be removed. Hansen 
& Lilje (1999) who used the original equation of the form (101) presented in Nagirner (1994) also noted the aforementioned 
misprints. 
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